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Abstract 

The response of a relativistic particle bound in a linear confining well is calculated as a function 
of the momentum and energy transfer, q, v. At large values of |q| the response exhibits scaling 
in the variable y = — |q|, which is proportional to the Nachtmann variable, ^. The approach to 
scaling is studied at smaller values of |q|. Scaling occurs at ~ |q| at relatively small |q|, and its 
validity extends over the entire ^ range as |q| increases; this behavior is observed in electron-proton 
scattering. About 10% of the response at large |q| is in the timelike region where v > |q|, and it 
is necessary to include it to fulfill the particle number sum rule. The Gross-Llewellyn Smith and 
Gottfried sum rules are discussed in the context of these results. 
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FIG. 1: The |q|-i^ plane. The spacehke region is above the |q| = v hne and timehke is below. 
Lines of constant > are parabolas which lie entirely in the spacelike region and approach 
q = as ^ cxD. The observed (Q^ > 0) response of the proton lies in the shaded area. 

Deep inelastic scattering (DIS) of leptons by hadrons is generally discussed in the frame- 
work of the naive parton model and the QCD-improved parton model using the operator 
product expansion, [m This approach has been very successful in determining the evolution 
of the structure functions as a function of the square of the four-momentum transferred to 
the hadron. In the leading order of the model the hadron is approximated by a collection 
of noninteracting quarks and gluons. The struck quark is assumed to be on the mass-shell 
both before and after its interaction with the electron. Based on this assumption, the lead- 
ing order response is predicted to be in the spacelike region for which the energy transfer v 
is less than the magnitude of momentum transfer, |q|, as a consequence of the inequality, 

V = ^|k + qP + m2 - y''|kp + m2 ^ |q|. (1) 

Here k and are the momentum and mass of the struck quark, respectively. The predicted 
response is discontinuous at the boundary |q| = z/ between space and timelike regions. 
The conventional variables of the parton model, = |qp — and the Bjorken x = 

2 



Q^/2M z/, used to describe the DIS structure functions of a hadron of mass M, are confined to 
the spacehke region of the |q|-z^ plane for positive values of accessible in lepton scattering 
experiments, as shown in Fig. |1|. Therefore we study the response, -R(q, v) as a function 
of V and |q| in the rest frame of the system [^], as is common practice in the many-body 
theory (MBT). Lines of constant |q| in Fig. |I] cross the photon line {v = |q|) and go into the 
timelike region. The observed {Q^ > 0) DIS response is limited to a narrow region in the 



|q|-z/ plane illustrated in Fig. |T]. It is bounded by the elastic limit, Uei = A/|qp + — M 
on one side, and by the photon line on the other. In the limit of large |q| the width of the 
observed response at fixed |q| is M. Lines of constant intersect the elastic limit curve at 
X = 1 and approach the photon line at small x. 

For a hypothetical scalar probe, the response is given by: 

i?(q, u) = 5^ K/| $^e^''-^ |0)p5(E, ~Eo-u) (2) 
/ j 

where is over all the particles and the over all energy eigenstates. It is viewed as 
the distribution of the strength of the state Ylj ^^^ ''^l^) o'^^^ ^^e energy eigenstates of the 
system having momentum q. It is not necessarily zero in the timelike, z/ > |q| region. 

The natural scaling variable in the MBT approach to DIS [0] is y = z/ — |q|. At large 
|q| the response is expected to depend only on y, and not on q and u independently. This 
variable is equivalent to the Nachtmann variable ^ since ^ 

e=^(|q|-^) = -^y. (3) 

In the limit of large the ^ = x, thus y scaling includes Bjorken scaling. However, both y 
and ^ span both spacelike and timelike regions at fixed |q| unlike x at fixed Q^. 

The particle number sum rule in MBT is obtained by integrating the response at large 
|q| over all z/ > 0: 

/ i j 
= 5]](0|e*'i-(''^-''')|0). (4) 
hi 

When q is large only the i = j terms in the above sum contribute, and therefore the integral 
gives the number of particles in the system. In contrast the sums of the response in the 
parton model are obtained by integrating the response over > at fixed Q^. These sums 




FIG. 2: ^ in GeV^ plotted against angular momentum £ for the first ten values of the radial 
quantum number n. The lines are linear fits to the calculated values. 

will fulfill the particle number sum rule only if the response in the timelike region is zero. As 
mentioned earlier, the response of a collection of noninteracting particles lies in the spacelike 
region. Interaction effects, however, can shift a part of the strength to the timelike region. 
Evidence for shifts caused by interactions is discussed in Ref. 

We have studied the exact response of a simple "toy" model which contains the basic 
features of relativity and confinement to obtain further insights on the possible response 
in the timelike region and it's effects on the sum rules. In this model we assume that the 
response of the hadron is due to a single light valence quark confined within the hadron by 
its interaction with an infinitely massive color charge. We model this interaction by a linear 
fiux-tube potential, and use the single particle Hamiltonian, 

H = ^\p\^ + ml + y^r (5) 

containing the relativistic kinetic energy operator. In the limit rriq = used here, the H can 
be cast in the form: 

i7 = a^/^(vW + r'), (6) 

where p' = p/a^^^, and r' = cr^/'^r are dimensionless. The response i?(|q|,z/) of the model 
then depends only on the dimensionless variables |q'| = |q|/cr^''^ and u' = vja^l^. The main 
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FIG. 3: The response for |q| = 10 GeV versus the variable, y = — |q|, for various Tq. 

conclusions of this work are independent of the assumed value of a; however, we show results 
in familiar units using the typical value ^/a^ = 1 GeV/fm. 

The model may be viewed as that of a meson with a heavy antiquark or that of a baryon 
with a heavy diquark. It is obviously too simple to address the observed response of hadrons. 
For example, it omits the sea quarks and radiative gluon effects contained in the DGLAP 
equations 0] to describe scaling violations. Nevertheless its exact solutions are interesting 
and useful to study scaling, the approach to scaling, and the contribution of the timelike 
region to sum rules. A similar model has been considered by Isgur et al. 0. 

The eigenstates of this Hamiltonian have the usual quantum numbers n, £, and m. They 
are expanded in the basis of j£(fc^r)Y^™(r). For each value of £ ^ 100 we consider 200 
(a = 1, . . . ,200) basis states each having je{k^^R) = 0, for R = 10 fm, and 300 basis states 
for R = 15 fm. The Hamiltonian is diagonalized within a sphere of radius R. The ground 
and relevant excited states have rather small radii, and are not influenced by this boundary 
condition. For example, the energies of the lowest 100 eigenstates of each i and the calculated 
responses at |q| ^ 10 GeV are essentially the same for R = 10 and 15 fm. Figure |^ shows 
a plot of the square of the energy eigenvalue for a given n versus i. The calculated energies 
lie on linear Regge trajectories having slopes within 10% of the classical estimate i^/a. 



5 



10 



^ 6 



o 



PS 



1 1 1 1 1 1 


1 Ir^iir, 1 

/ 


1 1 1 1 1 1 1 

-- q=3 GeV - 




/ ¥ 


q=4GeV 




/ \ 


-- q=5GeV _ 






- - - q=6GeV - 


- 




— q=10GeV " 


— ■ — 1 , , , 


> 

1 , , < 1 , , , 1 


, , , I'^r-T — , — 



000 



500 



1000 



-1500 -1000 -500 
y (MeV) 

FIG. 4: The response for values of |q| ^ 3 GeV versus the scahng variable, y = — |q|. 

The response is calculated including all states with n ^ 100 and i ^ 100 in Eq.(||). It is 
a sequence of 6 functions at z/ = E^^i — Eq. We verify that the full strength of the integrated 
response, 



(7) 



is obtained in this truncated basis for all values of the momentum transfer considered in this 
work with < 0.02 % error. In order to obtain a smooth response we assume decay widths 
for all the excited states dependent on the excitation energy u: 



r(z/ ^ Et) 



To (1 



-{u-Et)/Es 



(8) 



Here Et is the threshold excitation energy for meson emission below which the width is zero, 
and Eg parameterizes the approach to a constant width Fq at z/ > Et. We use Et = Eg = 100 
MeV and various values of Fq for illustration. Only the response at small values of |q| is 
sensitive to Et and Eg. The response including decay widths is given by: 

J(i|e'-'|0)r I ^^""^ ^ 
I ^ 



i?(q,z.) = ^|(/|e^'^-|0)p 



(9) 



27r J {Ei -Eo- uf + F2(z/)/4 ' 

The dependence of the response at |q| = 10 GeV on Fq is shown in Fig. |^. Since this response 
has v ^ Et and Eg, it depends only on Fq. 
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FIG. 5: The approach to scaling of the response for values of |q| ^ 2 GeV and |q| = 10 GeV 
versus the scaling variable, y = — |q|. 

Figure ^ shows the response calculated for values of |q| ^3 GeV as a function of y for 
To = 100 MeV. The scaling behavior is clearly exhibited; at large |q| the -R(|q|, v) becomes 
a function f{jj) alone. This scaling is equivalent to ^ scaling via (Eq. H). 

In Fig. ^ we show the response for Fq = 100 MeV, at various values of |q| ^2 GeV 
compared with that for |q| = 10 GeV, to study the approach to scaling. At small |q| the 
scattering is dominated by resonances, and the first inelastic peak is due to the lowest excited 
state with n = 1 and ^ = 1, 335 MeV above the ground state. In our toy model, the elastic 
scattering occurs at = or y = — |q|, since our hadron is heavy. This elastic scattering 
contribution is omitted from Fig. |. It has a strength of {0.85, 0.53, 0.26, 0.10, 0.010, 0.00086} 
for |q| = {0.25, 0.5, 0.75, 1.0, 1.5, 2.0} GeV. The energy dependence of F [Eq.(|D] implies that 
the inelastic response is zero for u < Et ot equivalently y < Et — \q\. 

For y ~ 0, i.e. for small ^, the response approximately scales at relatively small values 
of |q|, comparable to a^^^. As |q| increases, the range over which scaling occurs is extended 
to more negative values of y, i.e. to larger values of C,. The contribution of each resonance 
shifts to lower y and decreases in magnitude following the -R(|q| — > oo,y). This behavior is 
seen in the experimental data on the proton and deuteron and interpreted as evidence for 
quark-hadron duality. Thus the toy model seems to describe some of the observed properties 
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of the DIS response of nucleons. It exhibits y or equivalently ^ scahng at large |q| as observed 
0, and an approach to C, scahng similar to that seen in recent experiments. 

The -R(|q|, z^), and therefore the f{y) extend into the timelike {y > 0) region. The sum- 
rule given by Eq.(|^, counts the number of particles in the target. It is necessary to integrate 
over the timelike region to fulfill this sum rule. In the limit Fq = about 9.6% of the sum 
is in that region independent of ^/o^. It increases to 13.7% for Fq = 100 MeV and ^/o' = 1 
GeV/fm. The response expressed as R{Q^,^) also scales at large where |q| is necessarily 
large. It becomes a function of ^ alone. However, the integral: 

/■oo ;>|q| 

/ R{Q^ oo,^)d^= i?(|q| ^ cx),z/)rfz/ < 0.904, (10) 
Jo Jo 

because the contribution of the timelike region is omitted. The maximum value, 0.904 of the 

integral is obtained for Fq = 0. Here we have defined ^ = |q| — z/ without the conventional 

1/M scale [Eq.(0)]. 

The baryon number, or Gross-Llewellyn Smith (GLS) [|l], §| sum, Sgls is obtained by 
integrating the neutrino scattering responses over the range < a; < 1. In absence of 
any timelike response, in the limit oo, Sgls = 3, the number of quarks minus the 

number antiquarks in a nucleon. The GLS sum has gluon radiative corrections that have 
been calculated at next-to next-to leading order {i.e. up to 0{al)) using pQCD 0. The 
experimental values of Sgls obtained from revised CCFR |jTO[ xF^ are 2.55 ± 0.06(stat) ± 
O.lO(syst) and 2.80 ± 0.13(stat) ± 0.18(syst) at = 3.16 and 12.59 GeV^, respectively. 
They have been used to determine asiQ"^) assuming that the difference between 3 and the 
above experimental values is entirely due to gluon radiative effects. 

Recent measurements of the Drell-Yan cross section ratios, cr(p + d)/a{p + p) from E866 
at Fermilab |Tl|, [1^ determine the ratio of d{x)/u{x) and, in turn, the integral, 

f rfx[M(x) = -0.118 ±0.012 (11) 
Jo 

to a high level of precision. The valence isospin, or Gottfried, sum 

dr 12/"^ 
Sg= — [Fiix) - F,^ix)] = - + - dx [uix) - dix)] (12) 
Jo ^ J <J Jo 

has been measured by the NM Collaboration ||l3l. The analysis of their complete data set 

gives Sg = 0.216 ±0.027 [|14[ |l5l essentially independent of on the interval 0.5 < < 10 
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GeV^. Using these measurements we find a contribution due to valence quarks of, 

2 



Sg- 



[ dx[u{x) -d{x)] =0.294 ±0.030, (13) 
Jo 



about 10% below the expected result of |, with comparable error. 

The present work suggests that these sums over the spacelike region could be smaller 
than the theoretical expectation due to some of the response being shifted into the timelike 
region. Such a shift is due to the bound nature of the quarks in the hadron and is therefore 
a nonperturbative effect. 

The prevalent models of valence quark structure functions predict a large valence 
quark response in the spacelike region adjacent to the photon line which diverges as ~ x~^'^. 
Even though the x < 0.1 region is dominated by sea quarks, contributions from this region 
make up a large part of the experimental sums; about 50% of Sgls HH and about 30% of 



Sq fT^l- Response having ^ < 0.1 is within < O.IM ~ 94 MeV of the photon line. It is 



possible that there is some response beyond the photon line, moved into the timelike region 
due to nonperturbative effects associated with binding. 

It is known from e^e~ annihilation experiments that the vacuum has timelike response, 
and the nucleon presumably couples to it. In the toy model we have ignored this problem 
in order to focus on the possibility that the nucleon can have timelike response due to the 
bound quarks in the nucleon being off the mass-shell. 
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